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Abstract 

A new computational idea for continuum quantum Field theories is outlined. 
This approach is based on the lattice source Galerkin methods developed by 
Garcia, Guralnik and Lawson. The method has many promising features in- 
cluding treating fermions on a relatively symmetric footing with bosons. As a 
spinoff of the technology developed for "exact" solutions, the numerical methods 
used have a special case application to perturbation theory. We are in the pro- 
cess of developing an entirely numerical approach to evaluating graphs to high 
perturbative order. 

1 General Approach 

Conventional numerical evaluation of quantum field theory involves evaluating the path 
integral on a spacetime lattice using Monte Carlo integration methods. Monte Carlo 
methods have been very successful for many problems and have served to at least 
roughly confirm the standard model of particle interactions. However, many serious 
issues are either not yet solved by this approach or not likely to be calculable at all. 
These include dealing with actions that are not manifestly positive definite or which 
have important effects from the details of fermionic interactions beyond the quenched 
approximation. Consequently, it is worthwhile to search for supplementary methods of 

numerical calculation. 
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The source Galerkin method, although at a relatively primitive stage of devel- 
opment, appears to work on both the lattice and the continuum while dealing with 
fermions as easily as bosons. Of course, the continuum version of source Galerkin com- 
pletely avoids the usual lattice problem of fermionic multiplicity and the tricks used 
to minimize errors from this problem. Furthermore, our "source Galerkin" method is 
less restrictive as to the class of allowed actions. Source Galerkin tends to use signif- 
icantly less compute time than Monte Carlo methods but can consume significantly 
larger amounts of memory. The formulation of a source Galerkin approximation to a 
given field theory can also require much more initial work and input to the calculation 
than the usual Monte Carlo approaches. 

We confine this discussion to the continuum analysis of 4 interactions; examples 
of lattice calculations have been given elsewhere. JI], 0, |3[] Our continuum results to date 
include non-linear sigma models and four-fermion interactions. We have gauge theory 
calculations in progress. 

We start with the differential equations satisfied by the the vacuum function Z for 
a scalar field with interaction ^0 4 /4 coupled to a scalar source J(x) which satisfies 
the equation: 



The source Galerkin technique is designed to directly solve functional differential equa- 
tions of this type. Before we proceed to outline a solution technique, it is essential 



is a set of third order differential equations at every point in space time, it possesses, 
a priori, an infinite solution set. The solution studied in this talk will be the usual 
one which corresponds to the symmetry preserving solution obtained from evaluating 
a path integral with real axis definitions for the regions of integration. This is the 
solution that is regular in the coupling, g, as it approaches zero. 




(1) 



to point out that this equation by itself does not uniquely specify a theory. Since it 
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If we write (]T|) in the form: 

EjZ(J) = (2) 

The source Galerkin method is defined by picking an approximation Z*(J) to the 
solution Z (J) such that 

EjZ*(J) = R (3) 

where R is a residual dependent on J. We pick the parameters of our approximation to 
make this residual as small as possible on the average. To give this statement a mean- 
ing, we define an inner product over the domain of J: i.e. (A, B) = J dfi(J) A(J)B(J). 
In addition, we assume we have a collection of test functions which are members of a 
complete set: {(fi{J)}- The source Galerkin minimization of the residual R is imple- 
mented by setting the parameters of our test function Z*(J) so that projections of test 
functions against the residual vanish so that | \Z* — Z\ | 2 — > as the number of test 
functions — > oo. 

The equations defining the quantum field theory are differential equations in the 
field sources and spacetime. We can deal with the continuum by taking advantage of 
our knowledge of functional integration to evaluate integrals of the form: 

1 = J{dJ}e W [-J 2 (x)/e 2 ]P(J) (4) 

and use this to define an inner product of sources on the continuum as follows: 

(t( \ t( \ rt \ u \\ / en+m M^i---^i---?/m} n + meven 
{J{x l )---J{x n ),J{y 1 )---J(y m ))j=< . (5) 

I U otherwise 

where we have absorbed a factor of 2 by redefining e. 5 + is defined by 

5 + {xa(3 ■■■} = 5(x- a)5 + {(3 ■ ■ ■ } + 5(x - [3)5 + {a •••} + •••, (6) 

5 + {xa} = 5(x — a). (7) 
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In addition to this inner product definition, we need good guesses for approximate 
form for Z* and numerical tools to calculate, symbolically or numerically, various func- 
tions and their integrals, derivatives, and so on. For most of our calculations, we have 
found it very useful to choose a lesser known class of functions, with very suitable 
properties for numerical calculation, known as Sine functions. We take our notation 
for the Sine functions from Stenger:|| 



Sine approximations satisfy many identities, |J which make these functions very easy 
to use for Galerkin methods, collocation, integration by parts, and integral equations. 

With the definition of a norm and set of expansion functions, we can postulate an 
ansatz for Z: 



This ansatz with only free field structure is very simple and one would not expect 
to describe the entirety of a given theory with it. However, the strength of the Galerkin 
method to force convergence to correct answers is often strong enough to yield remark- 
ably good results for the lowest state of a theory. Extremely good approximations 
to theories are obtained as additional source structure combined with the appropriate 
spacetime structure is added to the ansatz. 

In principle, it is possible to allow the functions multiplying the sources to be totally 
arbitrary and then to fix them by expanding in Sine functions and apply the Galerkin 
technique. The associated computations soon demand more memory than is available 
from current computing systems. Instead we can use our knowledge of the spectral 
representations of field theory and graphical approaches (developed from perturbation 
theory but we emphasize that our approach is not perturbative) to introduce a beautiful 




tt(x — kh)/h 





(9) 



4 



and intuitive approach of regulated Lehmann representations to produce candidates for 
Z* . These build the appropriate spacetime Lorentz structure into our approximations 
and make the operational cost of our numerical approach independent of spacetime 
dimension. Since any exact two-point function can be represented as a sum over free 
two-point functions an appropriate representation of the two-point function greatly 
facilitates this approach. We choose as the basis of our numerical solutions a regulated 
Euclidean propagator structure: 

e ip-x-p 2 / A 2 



A(m; x) = (dp) 



p2 _|_ m 2 



(dp) I dse ip - x - p2/A2 - s(p2+m2) 





ds 



TV 



s + 1/A 2 



d/2 



(10) 

"4(H-1/A*) (11) 



(2tt) ,o 

This cutoff in the Galerkin test function as well as the definition of inner product 
assures that we never have to address the issue of a divergence during the course of a 
calculation, just as in the case of lattice calculation where the lattice spacing serves as 
the cutoff. 

This integral can be approximated using Sine methods 



A(m; x) 



h 



N 1 

y J- 



(27t)" ' e 

v I k=~N 



71 


d/2 


exp 


U + l/A 2 . 



X 



-z k m - 



±(z k + 1/A 2 ) 
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Here z k = e . For practical purposes, a sum over eighty terms with A = 10 gives 
accuracy out to twelve digits — more than adequate for most hardware floating-point 
representations. Thus we have a form for a two-point scalar Green function, regulated 
by the scale A 2 with constant computational cost regardless of spacetime dimension. 
We can take derivatives explicitly or by construction: 



<9 2 A(m; x) = m 2 A(m; x) — S(x) 



(13) 



where S(x) = e x2a2 /^ 
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Figure 1: One dimensional 4> 4 mass gap versus coupling (dashed line gives exact from 
Hioe and Montroll, 1975) 

From this representation, we can directly construct a fermion two-point function: 

S(m; x) — (7 • d — m) A(m; x) (14) 

These representations mean that free scalar and free fermion results are exact and 
immediate in any Galerkin evaluation of these trivial cases. 

2 Results: </> 4 

We itemize some results obtained using a regulated single propagator with parameters 
set by the Source Galerkin method. At lowest order, our ansatz for the generating 
functional is 

Z* = exp J ~j x G xy j y . (15) 

Results for this ansatz are given in Figure [I]. These results are strikingly accurate and 
can matched up essentially exactly with results of Monte Carlo calculations in two and 
higher dimensions. 



(3 




Figure 2: Additional connector-based ansatzen for the four-point function, H, in A0 4 . 
In the bottom row, we have two contact ansatzen on the left, followed by two mediated 
ansatzen. 

We can enhance these results by including additional 4 source terms in Z*. Some 
simple additional terms that we include with weights and masses to be calculated using 
the Source Galerkin technique are the terms of the forms given in Figure |2|. The effect 
of adding a fourth order term is shown in Figure ^. 

3 Perturbation Theory 

The representation that we have given for the free propagator fllCf ) provides the basis for 
an entirely numerical expansion of any Feynman graph. We have examined the possibil- 
ity of of generating all graphs of a theory through automatic functional differentiation 
of a vacuum amplitude and then numerically analyzing the resultant expressions using 
Sine expansions. || [7j This method generates and evaluates graphs at a speed which 
to our knowledge far exceeds any method used in the past. We believe that it will be 
possible to calculate to quite high orders in perturbation theory. 
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Figure 3: Comparison of four- if approximation with lowest order and exact answer, in 
one-dimension (A 2 = 70). 

4 Conclusions 

We have discussed a method for numerical calculations for field theories on the con- 
tinuum; this method being based on the original lattice source Galerkin technique. 
We presented a Lorentz-invariant regulated representation derived from Lehmann rep- 
resentations. This approach has the computational advantages of minimal memory uti- 
lization and parallelizable algorithms and also allows direct representation of fermionic 
Green functions. Finally, a number of useful peripheral calculations can be made using 
this approximate representation: one can calculate diagrams in a regulated pertur- 
bation theory, as well as calculating dimensionally regularized loops numerically. In 
general, this technique of evaluating field theories takes advantage of the symmetries of 
the Lorentz group; future work includes the extension of the method to more general 
internal groups, such as gauge groups or supersymmetry. 
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